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Abstract 



The general equations describing hydrostatic equilibrium are developed for 
Non-singular Gravity. A new type of astrophysical structure, a Super Dense 
Object (SDO) or "Dark Star", is shown to exist beyond Neutron star held 
strengths. These structures are intrinsically stable against gravitational col- 
lapse and represent the non-singular alternative to General Relativity's Black 
Holes. 
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I. INTRODUCTION 



Recently it was discovered that a generalisation of Einstein's theory of General Relativity 
(GR) yields an analog of the Schwarzschild solution which is everywhere non-singular [|TJ. 
This remarkable result was found for a hitherto neglected sector of Moffat's Non-symmetric 
Gravity Theory (NGT) 0. The neglected degrees of freedom were found to affect drastically 
the strong field behaviour of the gravitational field. The smallest, non- vanishing contribution 
from this sector removes event horizons and renders the curvature everywhere finite. This 
sector of NGT has been dubbed Non-Singular Gravity theory or NSG in order to distinguish 
this new physics from the singular truncation of NGT studied in the literature to date. 
The results in Ref. QXJ] have now been generalised to include the analog of the Reissner- 
Nordstrom solution. It was found that the electric field is everywhere non-singular and both 
the electromagnetic energy density and spacetime curvature are everywhere finite ||. 

These results immediately raise several questions: (1) What do static structures look 
like? (2) Can they undergo gravitational collapse? (3) What is the end result of such a 
collapse? This paper is primarily devoted to answering the first of these questions, although 
tentative answers are proposed for all three questions. It will be shown that two distinct 
classes of interior solution exist. The first class is very similar to the usual GR solutions, 
while the second class describes a new and unusual Super Dense Object (SDO). The GR 
type structures are unstable against gravitational collapse, although less so than their GR 
counterparts. The SDOs on the other hand are intrinsically stable since their binding energy 
tends to zero as they become more compact. Stable SDOs only exist with radii in a small 
range near their "Schwarzschild radius" R = 2M. SDOs with radii much less or much 
greater than this value are unstable against gravitational expansion. For these reasons, 
it is conjectured that SDOs represent the non-singular alternative to GR's Black Holes. 
Astrophysical phenomena such as Quasi-Stellar Objects (QSOs) and (Active) Galactic Nuclei 
(AGNs) may well be SDOs. 

The outline of this paper is as follows: First the field equations are presented. From these 
the equations describing hydrostatic equilibrium are derived. The equilibrium equations are 
then studied in a small region near the origin, where the existence of two geometrically 
distinct branches of solution is established. Each branch of solution is then studied in turn 
using a mixture of analytic approximations and numerical surveys. Physical bounds on the 
"skewness parameter" , which controls the departure from GR, are obtained by considering 
Neutron stars. 

II. FIELD EQUATIONS 

Einstein's General Theory is based on the assumption that the metric of spacetime is 
a symmetric tensor. By removing this assumption, Moffat was able to construct a non- 
symmetric theory of gravity where both the metric and connection are non-symmetric. 

The NGT Lagrangian with sources is given by Q 
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where R^iT) is the generalized Ricci tensor, W p is a Lagrange multiplier, T refers to the 
torsion- free connection , T pv is the energy-momentum tensor and S p is the conserved 
current which gives rise to NGT charge. Square brackets denote anti-symmetrisation and 
units are chosen so that G = c = 1 throughout. The field equations that follow from ([[J) are 



9nv,o ~ 9 P v^% - 9 w Kv = 0, (2) 

(V^99 [H ), U = ^V^9S^ (3) 

= ~W M + 8vr(T Ml/ - \g^T). (4) 

where A^ 7 = + -D^ 7 , and the tensor D satisfies 

4 , . 

+ = ~o 7r5P (^^P ~ 9pp9au + 9pu9[ap]) • (5) 



The generalised Bianchi identities, which result from the diffeomorphism invariance of NGT, 
give rise to the matter response equations 

-^=1 (g m {V=9~T^) u + g w (^T^) J + T>» {g mv + g pv>tx - g mp ) + \w [p , u] S» = . 

(6) 

Taking the matter to be described by a perfect fluid, the energy momentum tensor is found 
to be U 

Tpv = (p + V)UpU v - pg pv , (7) 

where u v is the fluid's four-velocity and p and p are the usual internal energy density and 
pressure, respectively. The Abelian gauge invariance of the Lagrangian under W p — > W p +X iP 
ensures the conservation of NGT current density. This is evidenced in (|3|) since 

MV^S^ = W=gg^)^ = . (8) 

The NGT charge contained in a sphere of radius r, in spherical coordinates, is given by 

l 2 (r) = J S°y/=g~drd6d<j) . (9) 

The static spherically symmetric interior case was studied for the truncated version of 
NGT by Savaria ||. In this work the pure divergence free parts of g^ were dropped. Even 
with this truncation, the analysis was more complicated than the usual GR case due to the 
inclusion of additional source terms arising from the NGT charge. In what follows, matter 
will be taken to be NGT charge neutral so that = and l 2 (r) = 0. This condition 
explicitly enforces the strong equivalence principle. 

The most general, static spherically symmetric metric for NGT was found by Papapetrou 
0] to be 
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(10) 



Taking l 2 (r) to be zero demands w{ 



0. The field equations are further simplified by 



making the allowed coordinate choice (3{r) = r . Both of these simplifications will be made 
in the derivation of the hydrostatic equilibrium equations. 



III. EQUILIBRIUM EQUATIONS 

The field equations can be most compactly expressed in terms of the following quantities: 
A=- log(r 4 + f 2 ) , B — arctan(r 2 //) • (11) 

The usual GR expressions are recovered when A = 21ogr and B = n/2. With w(r) = 0, 
the field equation (|j) is automatically satisfied. This is a generic feature of NSG, as the 
divergence-free part of g^*"* is the NSG sector of NGT. It was these divergence-free degrees 
of freedom which were dropped in previous work on NGT. The (tt), (rr), (88) and 
components of (f§) yield 



(log 7 )"-i(log 7 )'(log^^ + A' (log 7)' =87ra(p + 3j9) , (12) 



lita(p-p) , (13) 



-2A" + A' (log a)' - ((A 1 ) 2 + (B'f) - (log T )" + i (log 7 )' (log 

(fB'-r 2 A'\ ,fr 2 B'-fA'\ (fB'-r 2 A'\, N , 2 , . . . 

( r 2 B' + fA' \ ul ( fB'-r 2 A' \ f r 2 B' + fA' \ n , 
q+ { 2a ) ~ B { 2a ) + { 4a ) {hgai) =~^(P-P) ■ ^) 

The constant q comes from the Lagrange multiplier field = 3qcos8/2, and prime denotes 
derivatives with respect to r. One of these equations may be replaced by the r component 
of the matter response equation (^), which simplifies to read 

p' = --(p + p)(log7)'- (16) 

This relation is identical to its GR counterpart. The agreement with GR follows from NSG 
obeying the strong equivalence principle when the NGT charge is zero. 

The set of equations (pT2|)- (p~5|) are best studied as three useful combinations. The first 
combination is / times (|14]) plus r 2 times (|i~5|) , which yields 



The second combination is / times flT5|) minus r 2 times fll4|), which yields 



(18) 

The third useful combination is (pT3|) plus (|T2|) minus 4a/ (r 4 + / 2 ) times (|T8|) which gives 



r 



2 -fq 1 /(r 4 + / 2 )A'V (A') 2 - (B') 2 - 2A" 



These three equations, along with the conservation equation ([If]), provide the most useful, 
complete set of field equations describing hydrostatic equilibrium in NSG. 

Starting with the first of these combinations, it is straightforward to verify that / = — qr 2 
solves equation (|l~7D . With this form for /, all the field equations collapse to be identical 
to their GR counterparts. Unfortunately, this simple solution is incompatible with the 
external Wyman solution J^fJ] (the NSG analog of the Schwarzschild metric) since having 
q non-zero leads to a jump discontinuity in the curvature invariant R^R^ at the edge 
of the structure. This follows from the fact that W<f, = for the Wyman solution, while 
W<f, = 3gcos6!/2 in the interior. Consequently, the matching of these solutions at the edge 
of the structure demands q = 0. 

The complexity of the field equations makes it difficult to combine them into a simple, 
fundamental hydrostatic equilibrium equation. However, considerable insight into the nature 
of the interior solutions can be gained by considering power-series solutions about the origin. 
Taking the following expansions: 



P = Po + Pi r + P2r 2 + ■ ■ ■ 


(20) 


P = Po + Pir + P2r 2 + ... 


(21) 


7 = 7o + 7i r + l2r 2 + ... 


(22) 


a = a + air + a 2 r 2 + . . . 


(23) 


f = fo + fir + f 2 r 2 . . . 


(24) 



and substituting them into ( |H|]l9l) leads to two distinct branches of solution. The first is 
characterised by having a = 1 and f = 0, while the second branch has a = and / 7^ 0. 
The first branch is very similar to the usual interior solutions found in GR, while the second 
branch is unlike anything seen before. 

Physically, the first branch can be thought of as a matter fields perturbing Minkowski 
space, while the second branch looks like matter fields perturbing the Wyman metric. The 
existence of this second branch will be seen to be intimately related to the non-singular 
nature of the theory. The fact that there are two branches of solution does not lead to any 
lack of uniqueness in the description of a structure since the two branches describe objects 
of very different field strengths. For example, the Wymanian branch of solutions cannot 
describe the sun, or a white dwarf or even a neutron star. Instead, this branch describes 
super dense objects (SDOs), which typically have radii of R < 2M. Each branch will now 
be considered in turn. 



5 



IV. MINKOWSKIAN SOLUTIONS 



The Minkowkian-type solutions have expansions about the origin which are very similar 
to their GR counterparts: 



f = ar 3 



2 2 5 o \ 

-vrpo - gVrpo + g a ) r 



P = Po~ g7r(Po + Po)(Po + 3p )r 2 + ... 
( 4 

7 = 7o f 1 + -7r (po + 3p ) r 2 + . . . 



a 




where 



M(r) = 47r 



21 , 
8n 



a' 
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(25) 
(26) 
(27) 

(28) 
(29) 



The above quantity is related to the gravitational mass beneath radius r. The relation is not 
exact, however, as the matching at the edge of the structure leads to a complicated relation 
between M.(R) and the mass of the structure M. The skew field, /, increases the mass of 
a structure above that of a structure with a similar matter distribution, p(r), in GR. From 
([29|) it is clear the mass increase results from the energy density a 2 of the skew field. Since 
a 2 has the dimensions of energy density, it is convenient to parameterise this quantity as 
some fraction of the central density: 



epo 



(30) 



Since the contribution to the effective energy density from / increases from 21<r 2 /4tt with 
increasing r, it is no surprise that even a small value of e causes significant changes. Nu- 
merical studies reveal that large departures from the usual GR results occur once e reaches 
tcrit ~ 0.1(2M/R). The fact that the critical value for e scales with the dimensionless "grav- 
itational potential" 2M/R is in keeping with the general physical picture of NSG. Large 
departures from GR occur when the gravitational potential is large. This relationship is 
made more precise by considering the matching of the interior solutions to the exterior 
Wyman metric. 

A. Matching Conditions 

The boundary between the interior and exterior solutions occurs at a radius r = R 
defined by p(R) = 0. The exterior Wyman metric is described by 

7 = e\ (31) 



a 



M 2 {v') 2 e-»{l + s 2 ) 
(cosh(au) — cos(6z/)) 2 



(32) 
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2M 2 e I/ [sinh(az/) sm.(bu) + s(l — cosh(a^) cos(6z/))] 
(cosh(azz) — cos(5z/)) 2 



where 



a = V 2 ' 6 = V 2 ' (34) 

and is given implicitly by the relation: 

e^(cosh(az/) — cos(fa/)) 2 2 = cosh(azz) cos(6z/) — 1 + ssinh(a^) sin(6z/) . (35) 

The dimensionless quantity s arises as a constant when solving the vacuum field equations || . 
Unlike the constant of integration M, which is related to a Gaussian integral of a conserved 
charge, the constant s is not connected with any conserved quantity in the theory. Unlike 
M, there is nothing in the theory which could explain why s might vary from one body 
to another. For this reason, it is reasonable to expect that s is a fundamental constant of 
nature. This conjecture is supported in the weak field regime where s is seen to be a coupling 
constant setting the relative strength of the skew and symmetric fields. In what follows, s 
will be taken to be a universal constant of nature. 

Formally matching the solutions at r = R is not very enlightening, since 7 is only known 
as an implicit function of r. More can be learnt by matching onto a small M/r expansion 
of the Wyman metric, where the metric functions take the near-Schwarzschild form: 

2M s 2 M 5 As 2 M 6 , , 

7 = 1 - — + 15^ + T5^ + --- ' (36) 

/ 2M 2s 2 M 4 7s 2 M 5 87s 2 M 6 V 1 

a= { 1 - — + ^^ + ^^ + ^57^ + ---J ' (37) 

„ sM 2 2sM 3 6sM 4 

f - — + — + ^ + — < 38) 

These expansions are only valid for 2M/r < 1, which is not a serious restriction seeing as 
most structures have 2M/ R < < 1 . 
Matching a and / at r — R reveals: 

c 2 M 4 / 1 K 77T- \ 

M = Mo+ ii& ( 11Mo + 4 " TE (po + 3po)R3 + ' ' ') ' (39) 



a 



2 _ s 2 M 4 . 4M 4vr 



9.R 6 

where M is the usual GR expression for the mass 
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l + ^ + ^(3po-5p ) J R 2 + ...) , (40) 



M = M GR = 4vr p r 2 dr . (41) 

JO 
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Since the corrections to the GR expression for the mass of the structure are always positive, 
the gravitational mass of a body in NSG will exceed that of a body with a similar density 
distribution in GR. This increase in mass is primarily due to a decrease in the binding energy, 
which occurs as a result of the repulsive contributions to the gravitational force coming from 
the skew sector. If s is too large, the repulsive contributions from the skew sector can cause 
a structure to become unbound. More moderate values of s serve to stabilise a structure 
relative to its GR counterpart. The leading term in fl4"0|) can be used to express e in terms 
of s: 



(42) 




where p is the "average density" of the structure, 



P 



3M 
AttR 3 



(43) 



The above expression relating e and s indicates that the critical value of s for which large 
departures from GR begin to be seen will scale as s cr u ~ R/M. This indicates that the 
most stringent bounds on s will come from the study of Neutron stars, which is no surprise 
considering Neutron stars are the only observationally verified astrophysical objects which 
require a relativistic description. The numerical results will show that if s exceeds Soo ~ 20, 
Neutron stars become unbound. 

One unusual feature of the matching at r = R is that the matching of the metric functions 
does not ensure the matching of the gradients of the metric functions. In GR, any pressure 
and density profile which smoothly approaches zero at the edge of the structure leads, via 
the field equations, to a smooth matching of a' and 7' once a and 7 have been matched. 
This in turn ensures that all curvature invariants smoothly match at r = R, which is the 
physically important condition. In NSG the curvature invariants can smoothly match at 
the edge of the structure even though there are jump discontinuities in the gradients of the 
metric functions, so long as these jumps obey certain relations which can be derived from 
the field equations. For example, the jump in 7' is given by the relation 



4R 3 



(/') 



'R 4 - p 
R 4 + f 



if] 



8R 3 



7 



2/ 



(44) 



where square brackets denote the jump in the enclosed quantity at r = R. In the GR limit, 
/ = and the above equation yields [7' ] = 0. In NSG there may or may not be a jump in 7', 
depending on whether there is a jump in /'. For the Minkowskian branch of solutions /' is 
always positive in the interior while /' is always negative in the exterior (if r > M), leading 
to a jump in /'. Equation fl4*4| ) serves as a useful check on the accuracy of the numerical 
solutions, and was found to be satisfied within error tolerances. 



B. Binding Energy 

The internal energy of a body is defined in the usual way to be E = M — m^N, where 
m^v is the rest mass of the N th constituent and N is the conserved particle number for this 
constituent: 
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(45) 



and n(r) is the proper number density for this constituent. The total internal energy can 
then be broken up into its thermal and gravitational components as E = T + V. The thermal 
energy is given in terms of the proper internal material energy density, e(r) = p(r) — mjvra(r), 
by 



R 



4tt v /«(r)(r 4 + /(r) 2 ) e(r)dr , (46) 



while the gravitational potential energy V is given by 

rR 



V = M — 4tt v /a(r)(r 4 + f{r) 2 ) p{r)dr . (47) 



The gravitational binding energy, Q = —V, is always smaller in NSG than it would be for 
the same mass distribution in GR. To leading order, the binding energy is given by 

Q = QGR _ 15 2 3 

8 

= Q Gn — eM— 

4tt \ p) 

= " - — (¥)--• (48) 

which explains why a large value of e (or s) causes structures to become unbound. For 
polytrope equations of state: 

p = K p( 1+1 /™) , (49) 
the Newtonian expression for the binding energy becomes 

5-n i? V ; 



which can be combined with ( |48"D to give an expression for the value of s for which the 
structure becomes unbound: 



72 V (R 



5(5 -n) y VM 



(51) 



The largest mass White Dwarfs are well approximated by a n = 3 polytrope, and have 
M/R pa 4 x 10 -4 , which bounds s to be below = 6.7 x 10 3 . The major source of error 
in (|5lD comes from the approximate expression for the binding energy, which can be off by 
a factor of two for White Dwarfs. The second term in (f48"D is essentially exact for White 
Dwarfs and agrees with the numerical results within 0.1%. A plot of the dependence of Q/M 
on s is displayed in Fig. 1. for a White Dwarf with M/R = 2 x 10 -5 . 
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FIGURES 



FIG. 1. (Q/M x 10 5 ) as a function of the skewness constant s for a White Dwarf with central 
density po = 1.45 x 1CT 13 km~ 2 . The lower line, (a), is the exact numerical result while the upper 
line, (b), is the analytic result employing the Newtonian approximation for the binding energy. 

The largest mass Neutron stars can be approximated by a n = 3/2 polytrope, and have 
M/R PS 0.1, which gives the far more stringent bound on s of Soo = 20. This value of will 
be roughly 30% too large due to the approximate value used for the binding energy. Since 
Neutron stars are known to exist, and since they are well described by GR, a tighter bound 
on s can be obtained by requiring that the repulsive force causes at most a 10% reduction 
in the binding energy. This then gives a maximum allowed value of s cr it = 6. Limiting NSG 
to be a 1% correction over GR for Neutron stars gives s crit = 2, and so on. 



— 1/2 

FIG. 2. f(r) in the interior of a Neutron star of radius R = 0.8322 p c = 12.5 km when 
s = 8.9. The upper line, (a), is an analytic approximation, while the lower line, (b), is the exact 
numerical result. 

Even if s is large enough to cause a 20% reduction in the binding energy, numerical 
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evaluation shows that the density profile is almost identical to that in GR, with the difference 
in p(r) never exceeding 0.1% anywhere in the structure. These studies were done using the 
equation of state 

5 

(52) 

where p c = l/(727r) km -2 is a critical density. Above this density the equation of state goes 
over to a n = oo polytrope. Choosing a central density of p — 0.2p c and taking s = yields 
a Neutron star with mass M = 0.0782, radius R = 0.8325 and binding energy Q GR = 0.0074 
in units of p^ 1 ^ 2 - From (51), these values lead to the prediction that s = 9.7 will cause 




a 20% reduction in the binding energy. The numerical results showed a 20% change was 
caused by having s = 8.9, and gave rise to a Neutron star with mass M = 0.0798 and radius 
R = 0.8322. Since the graphs of a, 7, p and p are essentially the same for both GR and 
NSG, the only function worthy of display is the skew field /. A plot of / is displayed in 
Fig. 2. for the above mentioned Neutron star. The exact numerical result is shown, along 
with the power series approximation given in fl2~5|) . Even for relativistic situations, such as 
Neutron stars, the approximate solution for / gives remarkably good results. 

It is interesting to note that the strongest bounds on the coupling of matter to NGT 
charge came from Neutron star calculations, where the NGT charge caused an increase in 
the binding energy ||. The tightest bound corresponded to a 25% increase in the binding 
energy. If both sectors of NGT are taken into account, this bound is removed since a 
moderate value of s can cause a compensatory 25% decrease in the binding energy. This 
possibility is worthy of further investigation. 

If s is a fundamental, dimensionless constant of nature, like the fine structure constant, 
it is not unreasonable to expect s to lie in the range 0.001 < s < 1. For this range of 
values, NSG has essentially no effect on the description of any known astrophysical object, 
save Black Holes. Since both the Chandrasekar and Oppenheimer-Volkoff limits for the 
mass of White Dwarfs and Neutron stars rely on quasi-equilibrium reasoning, and since the 
equilibrium description of these structures is unchanged, the usual sequence of gravitational 
collapse will occur. The departure of NSG from GR comes as the collapse continues, not to 
a Black Hole, but to the SDOs to be described in the next section. 



V. SUPER DENSE OBJECTS 

The second branch of interior solution is characterised by having a rich geometrical 
structure near the origin. The metric near r = is of the form 

7 = 70 + 7 2 r 2 + O(r 4 ) , (53) 
a = a 2 r 2 + C(r 4 ) , (54) 
/ = fo + / 2 r 2 + (9(r 4 ) , (55) 

so that the proper volume scales as V p ~ r 2 , while two-surfaces of constant r have surface area 
S = Anr 2 . The curvature invariants are all constant near the origin, and are proportional 
to (S/Vp) 2 . For the vacuum Wyman metric the leading coefficients are given by 
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+ 0(s 2 ) 



(56) 



470 



2 (l + ^)) > 
4 - f + , 2 + 0( S 3 )) , 



(57) 



M 2 s 2 



/o 




(58) 



where a small s expansion has been employed. It should be noted that the limit / — > 
is achieved by taking M — > 0, not s — > 0. The Wymanian interior solutions possess the 
same geometry near the origin, although the coefficients in the expansion are affected by 
the presence of matter. Interestingly, the ratio 47o/«2 = s 2 M 2 remains almost completely 
unchanged in the presence of matter. This serves to highlight a major technical difficulty in 
obtaining numerical solutions to the field equations for this branch of solution. The standard 
method for obtaining numerical solutions relies on power series solutions near r = to 
initiate the integration out from the centre of the structure. This is usually no problem 
since the region near r = is essentially Minkowski space, and no a priori knowledge of 
global parameters, such as the structure's mass, is required. For the Wymanian branch of 
solutions this is not the case since the mass of the structure must be known before it can be 
determined! 



An obvious solution to this quandary is to break with convention and integrate from 
the edge of a structure in toward the centre. This can be successfully accomplished in GR 
by specifying a mass and radius and integrating to find the corresponding central density 
(although unphysical choices for the mass and radius will give rise to negative central den- 
sities). The success of this operation in GR is predicated on the fact that smooth matching 
of a and 7 automatically results in the smooth matching of all gradients. In NSG the jump 
in the gradients must also be specified. For the Minkowskian branch of solutions non-zero 
jumps could not be avoided, making this method difficult to implement. Thankfully, the 
Minkowskian branch can be treated using the conventional method. For the Wymanian 
branch of solutions /' has the same sign in the interior and exterior solutions, so it is not 



FIG. 3. The density profile p(r) for a typical "neutron SDO". 
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imperative that there be a jump in /'. In the absence of a proof that the magnitudes of /' 
should also match, it is difficult to proceed. While there is no guarantee that the derivatives 
match, there is no doubt that choosing all gradients to match is one acceptable possibility. 
This assumption will now be made in order to proceed. 

A wide numerical survey was made in order to cover a range of equations of state and 
values of s. While the details differ, all SDO solutions studied did fit a standard picture. The 
generic features include a decrease in binding energy for radii above or below a central value, 
Rfi, which typically lay in the range 1.1 M < < 2.1 M. The central density exhibited 
similar behaviour, but always peaked at a slightly smaller value of radius, R p < Rq, and 
dropped more rapidly for larger radii. SDOs with radii slightly larger than R p had density 
profiles which grew and then dropped, while SDOs with radii below R p had monotonically 
increasing density profiles. A comparison of binding energies at fixed po for SDOs with 
radii on either side of R p , revealed that the SDOs with radii less than R p were more tightly 
bound and thus energetically favoured. This is in agreement with the usual expectation 
that stable structures must have monotonically increasing density profiles. Of course, usual 
expectations should not be relied upon for these unusual objects. Varying the value of s did 
not change this overall picture, although it did increase Rq and R p somewhat, and tended 
to reduce peak binding energies (which could reach huge values for small values of s). 



FIG. 4. The skew function f(r) inside a typical "neutron SDO". 

The most important of the above properties is the intrinsic stability of all SDOs against 
gravitational collapse. If one tries to crush a SDO into a smaller radius its binding energy 
decreases, and like a cosmic rubber ball, the SDO will spring back to its equilibrium radius. 
This intrinsic stability allows stable SDOs to exist for any mass. SDOs represent a non- 
singular endpoint to gravitational collapse, and provide nature with a sane alternative to 
GR's Black Holes. While a SDO resides deep inside a gravitational well, and can have a 
surface redshift so large that it is essentially a "black star" , the gravitational gradients are 
mollified to the extent that matter pressure is able to support the structure. 

In order to illustrate these general features, the density profile and skew field are displayed 
for a "neutron SDO" . The equation of state is designed to roughly model a relativistic core 
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coated with a crust of non-relativistic neutrons: 

P=fl". '""I (59) 

Up 3 > p< i • 

Here units have been chosen where p c = 1. The density profile shown in Fig. 3. is for an 
SDO with s = 0.1, M=1,R= 1.298. This SDO had a central density of p = 9.63 x 10 6 , 
and a binding energy of Q = 26.02. The skew field / is displayed in Fig. 4., and illustrates 
the primary difference between the Wymanian and Minkowskian branch of solutions. The 
Wymanian branch is characterised by having a large central value of /, with f ~ AM 2 . 
The field undergoes a small initial increase before dropping sharply toward the edge of 
the structure. The Minkowskian branch, on the other hand, has / starting from zero and 
sharply increasing. Interestingly, the energetically disfavoured SDOs with radii above R p 
have skew functions which start near zero and increase to a peak value before decreasing. 
This suggest a series of quasi-equilibrium configurations may exist between the Minkowskian 
and Wymanian branch of solutions, although time dependent collapse calculations will have 



to be performed to verify this conjecture. Work is currently in progress |TD] to study the 
evolution of SDOs through gravitational collapse in NSG. 

The aforementioned neutron SDO has a surface redshift of z = 2.09 x 10 5 , which suggests 
the name Dark Star is appropriate for these objects. As s tends to zero, the surface redshifts 
associated with SDOs tend to infinity, making Dark Stars as black as Black Holes. 

The neutron SDOs with s = 0.1 become unbound for radii less than R « 1.264. The 
dependence of binding energy and central density on radii was studied in detail for the stiff 
equation of state: 

p — pi . (60) 

This involved finding the density profiles for over fifty different radii in order to produce 
plots of fl(R) and po(R). 



FIG. 5. The binding energy H as a function of radius R for a SDO with a stiff equation of state. 
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For most choices of s and equation of state fewer radii were studied, since it was sufficient 
to check that the general picture was the same. In Fig. 5. the binding energy is plotted 
against radius for this stiff equation of state with s = 0.1. A similar plot of the central 
density versus R for these SDOs is displayed in Fig. 6. For radii beyond R = 2.04 a new 
complication enters the picture since the central density is zero and there is a region of 
empty space near the origin. This would require an inner and outer matching to the vacuum 
Wyman solution, although this possibility was not pursued since such annular distributions 
become unbound. 



FIG. 6. The central density po as a function of radius R for a SDO with a stiff equation of state. 

The above description of SDOs is at best a rough sketch, based on a wide, but by no 
means exhaustive, numerical survey. Since little is known about the behaviour of matter for 
the kinds of densities being considered for SDOs, it is difficult to choose a sensible equation 
of state. On the positive side, the generic features described for SDOs do not seem to be 
greatly affected by the choice of equation of state. The possibility that the gradients of a , 7 
and / exhibit jumps at the edge of some SDOs remains to be investigated. It is possible 
that a proof can be found to show the Wymanian branch of solutions never exhibit jumps 
in metric gradients. No work has yet been done to study the feasibility of describing QSOs 
or AGNs in terms of SDOs, although the prospects are very good. In short, much work 
remains to be done to understand these fascinating objects. 



VI. DISCUSSION 

The non-singular gravitational theory described in this paper reproduces every prediction 
of GR for weak to moderate gravitational fields to any desired accuracy (taking s small 
enough). Any experimental prediction based on systems with 2M/R < 1 will fail to tell 
NSG and GR apart, unless s is larger than ~ 1. If s is larger than this value there is some 
hope that detailed studies of Neutron stars may tell the two theories apart. 

The only unique signal for any value of s would come for the super dense objects believed 
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to exist at the centre of most galaxies. If these galactic nuclei have masses above 1O 8 M , 
GR predicts that a star passing by will be torn apart by tidal forces once it is within the 
horizon of the massive Black Hole, so the star would vanish without a trace. If the galactic 
nuclei was a SDO on the other hand, the star's death would be visible [|ll] . Unfortunately, if 
s is too small, the radiation given off as the star is torn apart would be so heavily redshifted 
that no signal could be seen above background. 

In the absence of experimental tests to tell NSG and GR apart, there is only the fol- 
lowing aesthetic distinction: Any non-zero value of s leads to a non-singular endpoint to 
gravitational collapse - a Super Dense Object or Dark Star. In contrast, GR is always fated 
to have a singular endpoint to collapse - a Black Hole. Why cling to the concept of a Black 
Hole when the smallest imaginable change to GR banishes these objects from our picture of 
the universe? 

Infinities belong in mathematics, not nature. 
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